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criteria  which  are  independent  of  the  basic  difference  scheme,  and  are 
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In  this  paper  we  study  one-step  finite  difference  approximations 
of  scalar,  well-posed,  hyperbolic  problems  in  the  quarter  plane 
X > 0,  t > 0.  We  assume  that  the  basic  scheme  is  explicit  and  dissipative 
(hence  stable  in  the  usual  Von  Neumann  sense),  and  that  the  complimentary 
boundary  conditions  are  of  a general  two- level  type. 

Kreiss'  theory  in  [5]  provides  a sufficient  stability  criterion 
for  the  above  approximation.  Briefly  speaking,  stability  is  asstired  if 
the  operator  which  represents  the  overall  numerical  algarithm  has  no 
eigenvlues  on  the  boundary  or  in  the  exterior  of  the  unit  disc.  !Hiat  is, 
Kreiss'  criterion  is  (indirectly)  stated  in  terms  of  the  basic  scheme  as 
well  as  the  boundary  conditions.  Our  analysis  uses  Kreiss'  theory  and 
provides  stability  criteria  which  do  not  take  into  account  the  basic 
scheme,  but  instead,  are  given  solely  in  terms  of  the  boundary  conditions. 

Ibis  is  done  for  the  outflow  as  well  as  the  Inflow  cases. 

We  begin,  in  Section  1,  by  discussing  the  (somewhat  simpler)  outflow 
problem  together  with  the  set-up  which  allows  us  to  quote  Kreiss'  criterion. 
This  leads,  in  Section  2,  to  our  first  boundary- dependent  conditions. 

In  Section  3 we  use  the  new  results  to  derive  simple  criteria  for  the 
t common  case  where  the  boundary  conditions  are  translator^;  i.e.,  where 

I boundary  conditions  at  different  grid  points  are  determined  by  the  same 

' coefficients.  Among  other  results  we  show  that  if  the  boundary  oondi^inns  arp 

generated  by  an  explicit  stable  scheme,  then  overall  stability  is  assured. 
Indeed,  we  find  it  very  easy  to  apply  the  criteria  obtained  in  this  section 
to  popular  boundary  treatments  such  as  extrapolation,  explicit  Euler,  impli- 
cit Euler,  and  the  Box-Scheme.  All  those  examples  are  found  to  be  uncon- 
ditionally satisfactory  in  the  sense  that  when  augmenting  arbitrary  (stable) 


L. 


3 


dissipative  schemes,  they  always  maintain  stability. 

Finally,  in  Section  4,  we  slightly  modify  the  previous  criteria 
in  order  to  fit  the  inflow  problem.  As  an  exan^le  we  show  in  this  case 
that  homogeneous  summation  at  the  boundary  preserves  stability. 

In  a forthcoming  paper,  we  discuss  an  extension  of  the  above  theory 
to  nondlsslpative  and  multi-step  approximations.  The  main  key  to  these 
cases  is  given  in  the  theory  of  Gustafsson,  Kreiss  and  Sundstrom,  [3]. 

It  is  hoped  that  the  present  work  can  also  be  extended  to  include 
systems  of  equations  where  outflow  and  inflow  unknowns  Interact  at  the 
boundary.  Such  extensions  should  be  possible  at  least  for  certain 
families  of  difference  approximations. 


The  Outflow.  Problem  and.  I^elss * Criterion. 

2 

Consider  the  scalar,  hyperbolic,  L (0,x)  well-posed  initial 
value  problem 


1 


(1.1  ) au(x,t)/3t  = adu(x,t)/dx;  a > 0;  x>  0,  t > 0;  u(x,0)  = f(x)  . 

To  approximate  (l.l),  we  introduce  a mesh  size  ^>0,  At>0;a  j 

grid  function  = v(vAx,t),  v = 0,  + 1,  + 2,  ...  ; and  a consistent  j 

i 

finite  difference  scheme  i 


(1.2) 


Here, 

fixed 

values 


V (t  + At)  = Qv  (t)  , V = 1,  2,  3,  ...  , 


Q = I)  a,E^;  Ev,,  = 

J=-r  J ^ 


vM-1 


the  a.  are  constants  depending  on  the  coefficient  a and  on  the 
J 

ratio  X a At/Ax,  such  that  a ^ and  a do  not  vanish.  Initial 

•r  p 

are  given  by 


mmmmr 


m 


k. 


\(0)  = fy  » V = 1,  2,  3,  ...  . 

Throughout  the  entire  paper  two  assun^tlons  will  be  taken  for 
granted. 


ASSUMPTION  1.  The  difference  scheme  is  dissipative,  i.e.,  for 
some  5 > 0 and  natural  u,  the  as^llficatlon  factor 


satisfies 


Q(l)  = S a.e^^^  , -7r<  |<  TT  , 

J=-r  J 


IQ(0I  < 1 - , -Tr<i<Tr  . 


Clearly,  in  our  scalar  case,  dissipativity  implies  that  ft(|)  is  power 
bounded  (by  1),  which  is  equivalent  to  the  (strong)  Von  Nuemann  stability 
of  the  scheme  should  it  be  applied  to  the  pure  Cauchy  problem  for  - «><  x<  «. 


ASSUMPTION  2.  The  operator  Q is  not  right-sided,  i.e.,  we  have 
r > 0. 

Actually,  using  dissipativity,  it  was  shown  in  Corollary  1 of 
[2]  that  Q must  be  two-sided  if  it  is  to  be  consistent  with 
du/dt  = adu/dx  for  arbitrary  a ^ 0. 

The  fact  that  r > 0 implies  that  in  order  to  approximate  the 
differential  problem  in  the  quarter  plane  x > 0,  t > 0,  the  difference 
scheme  has  to  be  augmented  by  numerical  boundary  conditions.  More  pre- 
cisely, we  have  to  specify,  at  each  time  step,  boundary  values  at  the 
r grid  points  p = 0,  -1,  ...,  -r  + 1.  We  do  it  by  setting  boundary 
conditions  of  the  form 


I 


We  choose  to  write  (I.3)  in  its  present  form  because  it  allows 


us  the  following  two  options  which  we  use  later  in  various  exan^les 
(i)  The  boundary  conditions  are  of  a one-level  type;  that  is 


vanish,  and 


(ii)  The  conditions  are  two- leveled,  i.e.,  not  all  the  c 


vanish,  and 


In  any  event  it  is  understood  of  course,  that  the  required  boundary 


values  must  be  computed  in  the  specified  order,  namely,  n = 0,  -1, 


The  difference  approximation  is  completely  defined  now  by  the  stable 
scheme  (1.2)  together  with  the  boundary  conditions  (I.3),  and  we  raise 


the  question  of  overall  stability.  For  this  purpose  we  consider  the  space 


which  satisfy 


of  all  grid  functions  w 


Upon  defining  inner  product  and  norm  by 


/p  becomes  a Hilbert  space  which  is  a discrete  analogue  of  L (0,«>) 


6. 


Now,  the  difference  approximation  may  be  written  as 

v(t  + At)  = Gv(t)}  v(t)  , v(t  + At)  € ig  , 

where  G ; /g  ■^2  ® linear  bounded  operator,  defined  by  the  basic 

scheme  (1.2)  together  with  the  boundary  conditions  (I.3).  Thus,  the 
numerical  algorithm  is  embeded  in  and  is  said  to  be  stable  if  there 

exists  a constant  K,  independent  of  At,  such  that 

• l|v'(t)||  < k||v(0)||  for  all  t * nAt  , m = 0,  1,  2,  ...,  and  v(0)  e ig  . 


The  above  formulation  with  Option  1,  is  equivalent  to  the  scalar 
version  of  Kreiss'  problem  in  [5]*  In  fact,  it  is  not  hard  to  see  that 
the  various  resvilts  of  [5]  hold  also  when  the  boundary  conditions  are 
those  of  Option  2 (e.g.  [2]).  We  assume  therefore  that  the  reader  is 
familiar  with  the  details  of  [5],  and  conclude  this  section  by  quoting 
Kreiss'  stability  criterion. 


THEOREM  1 (Kreiss).  The  approximation  defined  by  (1.2),  (I.3) 


stable  if 


(i)  the  corresponding  operator  G has  no  eigenvalues  z with 
|zl  > 1,  z 1 ; 

(ii)  z = 1 is  not  a (possibly  generalized)  eigenvalue  of  G. 


iaiy-Dependent  Stat 


'riteria . 


Associated  with  the  difference  scheme  (1.2)  is  the  characteristic 


equation 


(2.1) 


S a.K^  - z = 0 , 


7. 


whose  r + p roots,  k^,  are  continuous  functions  of  z.  Kreiss 
has  used  dissipativity  and  consistency  to  show,  in  Lennna  2 and  in  part 
of  the  proof  of  Lemma  7 of  [5], the  following. 

LEMMA  1 (Kreiss).  (i)  For  z with  [z|>l.  z 1.  the  r+p  roots 
of  (2.1)  split ; p with  [k^I  > l and  r with  0 < |k^|  < i. 

(ii)  Suppose  that  the  coefficient  a in  (l.l)  is  positive 
(negative),  and  let  z with  [zl  > 1,  z 1,  tend  to  1.  Then,  the 
above  splitting  is  maintained,  except  for  a single  which  tends  to 
1 from  the  exterior  (interior)  of  the  unit  disc. 

Now  let  z with  [z[  > !•,  be  given.  If  z is  an  eigenvalue  of  the 
operator  G,  then  there  exists  a corresponding  (nontrivial)  eigenvector, 
g € ig,  such  that  Gg  = zg;  thus  by  the  definition  of  G,  g must  first 
satisfy  the  resolvent  equation 


(2.2) 


^ 3’  ••• 


Equation  (2.2)  is  an  ordinary  difference  equation  with  constant  coeff- 
icients; hence  its  most  general  solution  in  ig  is 

n “i-l 

where  = i«^(z)»  1 < i < n,  are  the  distinct  roots  of  (2.1)  which 
satisfy  0 < [k^I  <1,  each  with  multiplicity  m^;  and  the  j^(v) 
are  arbitrary  polynomials  in  v with  degree [P,  . (v)]  = k.  The  ct-  . 
are  parameters  to  be  determined,  and  by  Lemma  1,  their  precise  number  is 


In  particular  we  choose 


which  leads  to  a general  solution  of  (2.1)  6f  the  form 

"'A  "(k) -1'“  ■ . 


To  determine  the  ^ we  remember  that  being  an  eigenvector  of 
G,  g must  also  satisfy  the  boundary  relations 

(2.^)  (Gg)^  = zg^  , ‘ H = 0 , -r  + 1 . 

The  definition  of  G at  the  boundary  is  given  by  (I.3);  hence  (2.4) 
bec(»ies 


(2.5) 


S cf,i^g.._^^  = E cf°^  g,,^^  , ^x  = 0 


-HJ  -U+J  - ‘^u+j  > ^ = - ,...,-r  + l 

Thus,  substituting  (2.3)  in  (2.5),  we  finally  obtain 


m_g-l 


A £ k - » • 


h=0,  •••,  "r  +1  , 


which  constitutes  a linear  homogeneous  system  of  r equations  in  the  r 
unknowns  j^.  Clearly,  g is  an  eigenvector  of  G if  and  only  if 
not  all  the  in  (2.3)  vanish,  that  is,  (2.6)  has  a nontrivial 

solution. 


At  this  point  we  associate  with  the  boundary  conditions  (I.3)  a set 


of  rational  boundary- functions 


(2.7)  (zc<«  - ^ 

which  are  uni(iuely  determined  by  the  boundary 
Since 


— 0 ^ m • • f — r "t"  1 f 

coefficients  c^^\ 


ak  j=o 

then  system  (2.6)  may  be  written  as 


n (z,K  ) 

(2.8)  S S — y-  ^ 


i=l  k=0  ak 


k 


H = 0 , . . . , -r  + 1 


It  follows  that  the  coefficient  matrix  of  this  system,  which  we  denote  by 


(2.9a) 


D - D(z  j j m^,...,m^)  , 


is  of  the  form 


(2.9b)  D 

where  the  B( z , , m^ 


(2.9c)  B(z,K^,m^)  = 


We  recall  that 


— [B(  z , Kj^,m^) , . . . ,B(  z , K^,m^)  ] , 

),  1 < i < n,  are  r x m^  blocks  given  by 


1 

N 

1 

Rq(z,k) 

R_^(z,k) 

a 

R_j^(  z , k) 

m^-l 

R_j^(z»'') 

' as 

• 

• 

ait 

a 

a 

L .J 

J 

g is  an  eigenvector  of  G if  and  only  if  (2.8)  has 


a nontrivial  solution,  i.e.,  if  D is  singvilar.  This  gives  us 


10. 


LEMMA  2 . Let  z with  I z I >1  be  given,  and  let  with 

0 < fK^f  <1,  1 < i < n,  be  the  corresponding  distinct  roots  of  the 
characteristic  equation  (2.1),  each  with  multiplicity  m^.  Then  z 
is  an  eigenvalue  of  G if  and  only  if 

det  D(z  ; ° • 

Defining  a partition  of  r to  be  any  set  of  integers, 
which  satisfy 

l<in^<in2<  •••  + nig  + • • • + = r , 

we  are  now  ready  to  state  our  basic  boundary-dependent  criterion. 

THEOREM  2.  The  approximation  (1.2),  (I.3)  of  the  initial  value 
problem  (l.l)  is  stable  if  for  any  z with  |z|  > 1,  any  partition 
{m^J^l  2£  ^ and  any  distinct  values  ^ith  0 < Jk^J  < 1, 

det  D(z  ; k^,...,k^  ; m^,...,m^)  / 0 

Proof.  Take  an  arbitrary  z with  !zl  > 1 and  let  = k^(z), 

1 < i < n,  be  the  corresponding  distinct  roots  of  (2.1),  each  with 
0 < (k^(  <1  and  multiplicity  m^.  We  may  assume  that  the  m^  are 
ordered,  so  by  the  hypothesis, 

det  D( z f , . . . , j m^ , . . . ,m^)  ^ D . 

Hence,  by  Lemma  2,  z is  not  an  eigenvalue  of  G.  Since  z was 
arbitrary.  Theorem  1 asstires  stability,  and  the  proof  is  complete. 
Theorem  2 is  simplified  when  the  boundary  conditions  are  of  the 


one-level  type 


In  this  case  the  associated  functions  (2.7)  are 


V'^  = 2 c y ^ , U = 0,  -r  + 1 , 

^ J=0 


and  D of  (2.9)  is  given  by 


(2.11a)  D ^ D( ...  j ...  >111^^)  ^ [B( » • • • >®(  **ji’**^n^ ^ 


(2.11b) 


RqC'') 


Rn(0 


B(»t,,m,)  = 


d * 

’ dK  * 


’•••’  m,-l 
dtt  ^ 


, 1 < ■«  < n 


The  matrix  D of  (2.11)  no  longer  depends  on  z,  hence  Theorem  2 
becomes  z- independent,  and  we  have 

COROLLARY  1.  The  approximation  (1.2),  (2.10)  is  stable  if  for  any 
partition  1 ^ and  distinct  values  with  0 < |k^|  < 1, 


det  D(  * . . , K f ...  ^ ^ • 


ranslatorv  Boundary  Conditions* 


In  practice,  the  boundary  values  v^,  u=0,  ...,-r+l,  are 
usually  determined  by  a repeated  procedure;  that  is,  the  boundary 


12. 


conditions  are  of  the  translatory  form 

' “‘O’ * 1 • 

where  the  coefficients  Cj  are  independent  of  |i. 

The  rational  boundary- functions  associated  with  (3.I)  are 

R (z,k)  = S (zc^^^  - , U = 0,  ...,  -r  + 1 . 

^ j=0  ^ J 

■ Consequently 

R^(z,k)  = k^^q(z,k)  , n=0,  -r+1  , 

so  the  r X r matrix 

D 3 D(  z j J ™ [®( z > > • • • jB(  z j ] 

of  (2.9),  is  given  by  the  r x m^  blocks 


Kq  ( Z , It ) 

Ro(z,k) 

k"^Rq(z,k) 

t“\(z,K) 

.,-1 

k~^Rq(z,k) 

• 

• 

a 

• 

• 

• 

bK  ^ 

k‘^\(z,k) 

''■^\(Z,K) 

»c"’'‘^\(z,k)_ 

1 < i < q . 

Ihe  fact  that  D is  determined  now  by  the  single  boundary-function 
Rq(z,k),  in^jlies  the  following  significant  simplification  of  Theorem  2. 

THEOREM  3.  The  approximation  (1.2),  (3.I)  of  (1.1)  is  stable  if  for 
any  z with  [zl  >1  and  k with  0<  Ik|  <1, 


13- 


(3.2)  RoCz^it)  3 E (zc^.^^  - ^ 0 . 

d=o  ^ ^ 

Proof.  Take  an  arbitrary  z with  UI  > 1,  a partition  ^ 

of  r,  and  distinct  values  k^,  1 < i < q,  with  0 < [k^I  <1.  In 
order  to  prove  stability,  it  suffices,  by  Theorem  2,  to  show  that 


(3.3)  det  D(  z f j • ,m^)  ^ 0 . 

For  this  purpose,  let 


(3.^) 


0 

L <r 


kJRq(z,k^) 


= 0 , 


be  a vanishing  linear  combination  of  the  rows  of  D (where  we  denote  the 
transpose  of  a vector  by  a prime).  The  vector  relation  in  (3.^)  consists 
of  r scalar  equations. 


0 

2: 

^=-r+l 


3"? 


1 < i < q 


, 0 < k < m^  - 1 


which  we  write  as 

(3-5)  ^ vr*^] =0  - 1 5 ' < . , 0 <.<«,- 1 

Since  0 < |k^|  <1,  then  by  hypothesis,  the  left  member  in  the  above 


brackets  satisfies 


ik. 


, l<i<q  . 


Thus,  expending  by  Leibniz'  rule  and  using  induction  on  k > 0,  we 
find  that  the  right  member  in  (3*^)  has  vanishing  derivatives  at 
K = K . , i ,e  • , 


d^  r y 

— r 2j  O'  k 

dK , L h=-r+l  ^ J , 


= 0 , l<-*<q  > 


We  conclude  that  the  polynomial 


P(it)  3 S a K 
• p=-r+l  ^ 


r+p-1 


which  is  of  degree  r - 1 at  most,  has  r roots;  k^,  1 < ^ < q, 
each  with  multiplicity  m^.  Hence,  P(k)  30  and  the  coefficients 
cr^  must  vanish.  Py  (3«^)»  therefore,  the  rows  of  D are  linearly 
Independent,  so  (3*3)  follows,  and  stability  is  established. 

As  was  realized  in  the  previous  section,  if  (3.I)  is  reduced  to 
the  one-level  case 


then  the  boundary-functions 


R (k)  = 2 c.K^  ^ = 0,  ...,  -r  + 1 , 

^ d=o  J 


cease  to  depend  on  z,  and  Theorem  3 provides  us  with 


COROLIiARy  2.  The  approximation  (1.2),  (3.6)  is  stable  if 
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IL(k)  a S c.K^  / 0 , V K with  0 < |k|  < 1 . 

J=0  J 

Before  tvimlng  to  examples,  we  remark  that  the  case  r = 1 
(which  is  important  if  only  because  it  applies  to  the  dissipative  Lax- 
Wendroff  scheme,  [6])  is  translatory  by  definition.  In  this  case,  however, 
the  results  of  Theorems  3 and  Corollary  2 are  given,  respectively,  by 
Theorems  2 and  Corollary  1 since  for  r = 1 the  matrix  D is  reduced 
to  the  scalar  function  R^. 


EXAMPLE  1.  Let  the  boundary  values  be  determined  by  an  interpolation 
polynomial  of  arbitrary  degree  s - 1.  More  specifically,  we  extrapo- 
late at  each  time  step  from  v-(t),  ...,  v (t)  to  v-(t);  then  from 

*••»  etc.,  to  obtain  the  one-level  transla- 

tory procedure 

(3*7)  v^(t)  = D > ^1  = 0,  ...,  -r  + 1 . 

1 

The  associated  boundary- function  is 

RoC**)  = 1 - E (®)(-i)^'"V  = (1  - K)®  j 

j=l  ^ 

thus  Rq(i')  ^ 0 for  ^ with  0 < |k|  < 1,  and  Corollary  2 implies 
the  stability  of  (1.2),  (3.7). 

This  result  was  stated  by  Kreiss  in  [4]  and  was  proved  later  by 
CJoldberg,  [l] . 

In  the  remaining  examples  of  this  section  we  analyse  two-level 
boundary  conditions.  Exan^les  2,  3 and  3 were  studied  by  Gustafsson  . 
et  al.  [3]*  find  Example  4 by  Skdllermo  [7].  In  the  above  references. 
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however,  the  authors  did  not  consider  general  basic  schemes,  and 
attention  was  restricted  to  specific  approximations  such  as  Lax-Wendroff , 
Leap-Frog  and  Crank-Nicolson. 


EXAMPLE  2.  Consider  boundary  values  which  are  defined  by  extrapolating 
along  the  mesh-digonal.  The  two-level  version  of  this  procedure  is 
of  zero-order  accuracy,  and  is  given  by 

(3*8)  v^(t  + At)  = v^^^(t)  , ^ = 0,  ...,  -r  + 1 . 

The  corresponding  boundary- fimction  is  simply 

Rq(z,k)  = z - k . 

Ifence , 

y 0 , V IzI  > 1 , Ik|  < 1 , 

and  by  Theorem  3 the  approximation  (1.2),  (3.8)  is  stable. 

Note  that  (3*8)  is  inconsistent  with  the  differential  equation  (1.1), 
unless  \ 3 At/Ax  satisfies  \a  = 1,  where  (3*6)  coincides  with  our 
next  exaiiq>le. 

EXAMPLE  3.  Here,  the  boundary  conditions  are  generated  by  the  right-sided, 
first-order  accurate,  explicit  Euler  scheme,  which  yields 

(3*9)  v^(t  + At)  = v^(t)  + Xa[v^+3^(t)  - v^(t)]  , \ = At/Ax  , h = 0,  ...,  -r  + 1 . 
We  have 

Rq(z,k)  = z - 1 + \a(l  - k)  , 


3 


I 
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where  In  order  to  secure  stability  of  the  basic  scheme  (1.2),  one  may 
assume  that 

0 < \a  < 1 . 

It  follows  that  for  all  z,  k with  Iz|  > 1,  |kI  < l, 

IRq(z,k)I  = [z  - \aK  - (1  - Xa)I  > \z\  - IXanI  - [l  - \al  = 

= |zl  - XalK|  - 1 + Xa  = ([zj-  1)  + Xa(l  - |k|)  > o . 

Thus,  the  hypothesis  of  Theorem  3 is  fulfilled,  and  stability  of  (1.2), 

(3.9)  follows. 

EXi\MPLE  4.  Take 

(3.10)  v^(t  + At)  - Xa[v^_j^^(t  + At)  - v^(t  + At)]  = v^(t)  , p = 0,  ...,  -r  + 1 , 
which  is  determined  by  the  right-sided  implicit  Euler  scheme.  Since 


Rp(z,K)  = z(l  + \a  - XaK)  - 1 , (Xa  > 0)  , 


then  for  Izj  > 1,  [kI  < 1, 

IBq(z,ic)|  > IzI  • jl  + Xa  - XaK|  - 1 > |l  + Xa  - XaK|  - 1 > 

> 1 + Xa  - |XaK|  - 1 > Xa(l  - |k|)  > 0 , 

* and  Theorem  3 in^lies  stability. 

Although  the  above  examples  follow  without  difficulty,  it  is  clear 
that  the  more  con^lex  the  boundary  conditions  are, the  harder  it  is  to 
verify  the  hypotheses  of  Theorem  3 or  Corollary  2.  This  is  particularly 
true  in  the  two- level  case 
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<3-^5  ^ ^ « - 0 * 1 . 

where  the  relevant  stability  criterion  in  Theorem  3 involves  two  indepen- 
dent variables,  z and  k. 

In  order  to  provide  a z- independent  alternative  to  Theorem  3> 
consider  the  boundary- scheme 


(3-12) 


S^^^Vy(t  + At)  = S^®^Vy(t)  , V = 0,  + 1,  ± 2, 

g(i)  3 s^^^(E)  = E c^^^E^  , i = 1,  2 , 


d=o 


i 


which  generates  (3.II)  by  restricting  v to  the  values  0,  -1,  -r  + 1. 

We  prove 

THEOREM  4.  Let  the  boundary- scheme  (3-12) 

(1)  satisfy  the  solvability-condition 


(3-13) 


S^^)(k)  a E ^0  , Vl^I  <1  ; 


d=0 


(ii)  be  'stable  (in  the  usual  Von  Neumann  sense.) 

Then  (1.2),  (3. 11)  is  a stable  approximation  to  (1.1). 

Proof.  Hie  anqplif icat ion  factor  of  (3*12)  is  given  by 

s(i)  = . Ul  < TT  . 


where 


§(^>(5)  = » -e  = 1,  2 . 


Since  (3.13)  implies 
(3-l'*a) 


1 


I ; ' 

f 

I - ; 

\ i 
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then  this  amplification  factor  is  well  defined,  and  by  hypothesis  (ii) 
it  satisfies 

|s:Ol  <1  , VI  , 

which  we  write  as 

(3.1^to)  |S^°^(|)|  < IS^^\|)|  . 

The  boundary- function  associated  with  (3. 11)  is 

(3-15)  IL(z,it)  = h = zS^^^K)  - S^°\k)  , 

" J=0  ^ J 

So,  for  |z|  > 1,  we  use  (3.1^)  to  find  that 

|li„(z,e‘®)|  = |zS<''>(e‘b  - = |zS<^^0  - §<‘’^5)1  > 

>1^1  • |s‘^^(s)l  - |s(‘’\e)l  >0  . 

That  is,  the  equation 

(3.16)  Rq(z,k)  = 0 with  Iz]  > 1 , 

has  no  roots  k with  |k|  = i.  Since  the  roots  of  (3.I6)  are  conti'uous 
functions  of  z,  we  conclude  that  the  number  of  k with  jal  < 1 is  fixed 
for  (z|  >1,  and  can  be  determined  by  considering  large  values  of  |z|. 
Writing  (3*16)  in  the  form 

S^^^(k)  - z"^S^°\k)  = 0 , 

We  let  |z|  -»<»  and  use  (3.I3)  to  find  that  (3.I6)  has  no  roots  in  the  unit 
disc.  In  other  words,  if  [zf  > 1 and  Rq(z,k)  = 0,  then  a must  satisfy 
Ik)  >1.  By  continuity  therefore,  if  )z|  > 1 and  Rq(z,k)  = 0,  then 
|k(  ^1}  i.e., 

^ ^ , V Iz|  > 1 , Ik)  < 1 . 


J 


This  implies  (3*2)>  and  Theorem  3 completes  the  proof. 

We  note  that  the  stability  properties  of  the  boundary- scheme  are 
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often  known  In  advance , so  in  applying  Theorem  4 it  remains  to  verify 
only  the  solvability  condition  (3.I3).  Moreover,  if  the  boundary- 
scheme  is  explicit,  i.e.,  the  boundary  condition  are  of  the  form 


(3*17)  V (t  + At)  = L c V (t)  , n = 0,  ...,  -r  + 1 , 

^ J=0  ^ ^ ^ 

then  (3*13)  is  fulfilled  automatically,  and  Theorem  4 amounts  to 


CORQLLARf  3*  (3 ‘IT)  is  generated  by  a stable  scheme  (in  the 

usual  Von  Neumann  sense),  then  (1.2),  (3. 17)  is  stable. 

Reviewing  Exanqple  2 in  light  of  our  new  results,  it  is  clear  that 
the  generating  boundary- scheme. 


Vy(t  + At)  = v^j^(t)  , V = 0,  + 1,  ± 2,  ..., 


is  unconditionally  stable,  hence  Corollary  3 assures  stability. 

In  Example  3>  'tlie  explicit  Euler  scheme  is  known  to  be  stable  for 
Xa  < 1,  and  Corollary  3 again  implies  stability. 

Finally,  the  implicit  Euler  scheme  of  Exaii5)le  4 is  unconditionally 
stable,  so  to  comply  with  Theorem  k,  one  has  to  verify  that 

S^^^k)  a 1 - Xa(K  - 1)  ^ 0 , V Ik|  < 1 . 

Since  Xa  > 0,  then  k with  |k|  < 1 gives 


Re  s''  ''(k)  = 1 + Xa(l  - Re  k)  > i , 


and  the  result  holds. 

In  concluding  this  section  we  present  a somewhat  more  Intricate 
exaiiq;>le  due  to  Gustafsson  et  al.  [3]. 

EXAMPLE  3.  We  use  the  Box-Scheme  to  generate 


The  first  deviation  from  the  previous  discussion,  occurs  in  the  argu- 
ment leading  to  Lemma  2:  The  case  |zl  > 1,  z ^ 1,  goes  over  unchanged; 
hut  as  z -♦  1,  Lemma  1 implies  that  one  of  inner  roots,  of  the 

characteristic  equation  (2.1),  tends  to  1 as  well.  According  to  Kreiss* 
theory  this  root  has  to  be  incorporated  in  the  solution  g of  (2.3),  thus 
g may  fail  to  belong  to  ig,  and  z = 1 is  then  a generalized  eigen- 
value of  the  operator  G.  As  a resvilt.  Lemma  2 takes  the  form 

LEMMA  2'.  (i)  Let  z with  jz]  >1,  z ^ 1,  be  given,  and  let 
with  0 < |k^I  < 1,  1 < i < n,  be  the  corresponding  distinct  roots 

of  (2.1),  each  with  multiplicity  m^.  Then  z is  an  eigenvalue  of  G 
if  and  only  if  the  matrix  D of  (2.9)  satisfies 


det  D(z  ; k^,...,k^  ; m^,...,m^)  = 0 . 

(ii)  Let  K.  = K„(l),  2<  £ < n,  be  the  roots  of  (2.1)  which  corresponds 

£1  Ml  “ ~ . 

to  z = 1,  each  with  multiplicity  m^.  Then,  z = 1 is  a generalized 
eigenvalue  of  G if  and  only  of 

det(l  ; 1, J 5 =0  . 

Having  Lemma  2' , the  inflow  version  of  Theorems  2 and  3 follow 
easily  along  the  line  of  the  original  proofs.  We  first  state 

THEOREM  2'.  The  approximation  (1.2),  (1.2)  of  the  initial  boundary  value 
problem  (^.l)  is  stable  if  (i)  for  any  z with  Izf  >1,  z 1, 
aav  partition  2£  and  arbitrary  distinct  values 


with  0 < kJ  <1, 


det  D(z  ; k^,...,k^  ; m^,...,m^)  ^ 0 ; 
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(ii)  for  any  partition  of  r - 1 and  distinct 

with  0 < |k^|  <1, 

det  D(1  ; 1,k^,  . . , ,k^  ; 1, nig,.  ..,111^)  / 0 

As  before,  if  the  boundary  conditions  are  of  the  one-level  type 
(2.10),  then  D is  given  by  (2.11),  and  we  get  a z- independent  version 
of  Corollary  1 which  is  left  to  the  reader. 

Next,  in  the  translatory  case,  we  replace  Theorem  3 by 

THEOREM  3'.  The  approximation  (1.2),  (3.I)  of  (4.1)  Is  stable  if 
the  boundary-function  in  ( 3.2)- satisfies  (i)  Rq(z,k)  ^ 0 for  all 
fz|  > 1,  0 < |Kf  < 1;  (ii)  Rq(1,1)  ^ 0.  . 

REMARKS,  (a)  If  r = 1,  then  Lemma  1 suggests  that  k = 1 is  the 
only  relevant  root  for  z = 1.  Hence,  hypothesis  (i)  of  Theorem  3'  can 
be  restricted  to  z with  |z|  > 1,  z ^ 1. 

(b)  Had  the  boundary- scheme  which  generates  (3.I)  been  consitent  with 
(4.1a),  we  would  have  had  the  zero-moment  condition 

s . s =(°)  , 

j=0  J j=0  J 
which  is  equivalently  written  as 

Ro(lA)  = 0 . 

We  see  that  if  hypothesis  (ii)  of  Theorem  3'  holds,  then  the  boundary- 
scheme  must  be  inconsistent,  and  in  fact  may  have  no  acctiracy  with  respect 


to  the  differential  equation. 


2k. 


I 

I 

i 

f 

i 


1 

i 

I 

i 

! 


(c)  If  (3.1)  reduces  to  the  one-level  case  (3.6),  then  Theorem  3'  sim- 
plifies and  we  obtain 

COROLLARY  2'.  The  approximation  (1.2),  (3.6)  is  stable  if  Rq(i^)  ^ 0 
for  all  K with  0 < (k(  < l or  k = 1. 

EXAMPLE  6.  Let  the  boundary  conditions  be  of  the  simple  form 

Riese  homogeneous  sums  are  suitable  in  case  that  the  analytic  conditions 
(U.lb)  are  given,  for  instance,-  by  u(0,t)  =0  or  3u(0,t)/dx  = 0.  We 
have 

R (K)  = h I 

d=0 

so  RqC*')  / 0 for  all  ^ 0,  and  again,  by  Corollary  2',  the 
stability  of  (1.2),  (^.2),  is  assured.  This  last  result  is  also  due  to 
Kreiss,  [4]. 

Last,  we  consider  the  boundary  conditions  in  (3. 11),  and  state 
the  following  inflow- counterpart  of  Theorem  4. 

THEOREM  4',  Let  the  boundary- scheme  (3. 12) 

(i)  stalsfy  the  sol-yabillty-criterion  S^^\k)  ^ 0 for  all  (k|  < 1; 

, (ii)  satisfy  the  inconsistency-condition  Rq(1,1)  / 0; 

(iii)  be  stable (in  the  usual  Von  Neumann  sense.) 

Then  (1.2),  (3. 11)  is  a stable  approximation  of  (4.1). 

Proof.  Under  hypotheses  (i)  and  (iii),  the  proof  of  Theorem  4 implies 


that 

(^.3) 


^Qiz,K)  / 0 , V |zl  > 1 , < 1 . 


Hypothesis  (ii)  (which  in  paractice  can  be  easily  checked  by  inspection) 
together  with  (4.3),  constitute  the  assumptions  of  Theorem  3',  and 
stability  follows. 

We  conclude  by  recalling  that  when  the  boundary  scheme  is  explicit 
and  the  boundary  conditions  are  given  by  (3«17)>  hypothesis  (i)  of 
Theorem  4'  becomes  redundant. 
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